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Abstract. We have studied the influence of far ultraviolet (UV) radiation (6< hv < 13.6 eV) from a massive young stellar object 
(YSO) on the chemistry of its own envelope by extending the models of Doty et al. (2002) to include a central source of UV 
radiation. The models are applied to the massive star- forming region AFGL 2591 for different inner UV field strengths. Depth- 
dependent abundance profiles for several molecules are presented and discussed. We predict enhanced column densities for 
more than 30 species, especially radicals and ions. Comparison between observations and models is improved with a moderate 
^ . UV field incident on the inner envelope, corresponding to an enhancement factor Go ~ 10-100 at 200 AU from the star with 

f^ ' an optical depth t » 15-17. The chemical networks of various species are explored. Subtle differences are found compared 

^^ , with traditional models of Photon Dominated Regions (PDRs) because of the higher temperatures and higher gas-phase H2O 

10 ' abundance caused by evaporation of ices in the inner region. In particular, the CN/HCN ratio is not a sensitive tracer of the 

^^ ' inner UV field, in contrast with the situation for normal PDRs: for low UV fields, the extra CN reacts with H2 in the inner dense 

''^ , and warm region and produces more HCN. It is found that the CH+ abundance is strongly enhanced and grows steadily with 

JN ' increasing UV field. In addition, the ratio CH^/CH is increased by a factor of 10^-10^ depending on the inner UV flux. High-/ 

lines of molecules like CN and HCN are most sensitive to the inner dense region where UV radiation plays a role. Thus, even 
r^ . though the total column density affected by UV photons is small, comparison of high-y and low- J lines can selectively trace 

CL(' and distinguish the inner UV field from the outer one. In addition, future Herschel-HIFI observations of hydrides can sensitively 

JL I probe the inner UV field. 

H 
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1 . Introduction and drive a rich high temperature chemistry. However, there are 

/\ ■ also massive YSOs which do not show this molecular complex- 

Lj ' A newly formed massive star in the earliest stages of evolu- .. u <. u »• u »n u • 1 

!3 . -' ° ity, but where energetic phenomena can still cause chemical 

^ tion is surrounded by a dense envelope and is deeply embed- , tjju* u • j- ..c.u 

■ - - ■' ^ ^ ■' changes. Indeed, chemistry may be a unique diagnostic of the 

ded in its natal molecular cloud. Although highly obscured, 1 .• • ^ ^ j-a: . i, • 1 tu * 

° ° -' relative importance of different physical processes that occur 

these objects show signs of different stellar activity. Radio con- .,,. ,, -i/c/-. 1 j j j . j- c u 

■I '=' ■> within the YSO envelope, and a good understanding of chem- 

tinuum observations reveal compact Hii regions close to the ,. ,.1, ^ • 1 -. u a tu i,- a, a 1 »• 

^ ° istry is therefore crucial to comprehend the birth and evolution 



young star (e.g., Churchwell 2002) and molecular outflows n 
-' ° ^ ° ' ' of massive stars 

have been observed toward numerous massive star-forming re- 
gions (Shepherd & Churchwell 1996; Beuther et al. 2002a). 
Outflows and high energy photons from the YSO can strongly 
affect both the physical and chemical structure of the enve 



Much effort has been put into the quest of understanding 
the physical and chemical structure of massive YSO envelopes 
in recent years. Van der Tak et al. (1999, 2000), Mueller et 



, „, ^- 1 J .c • i^n,^ • al. (2002) and Beuther et al. (2002b) investigated the physi- 

lope. Observational studies of massive YSO environments re- , . , , , , . 

, . , , , . -.u 1 ■ . c ■ t cal structure of the envelopes around several massive young 

veal rich molecular spectra with a large variety of species (see ,^„ ,^s «tt , ■ ,■ ,- , ^r^ , , 

„, , , , .„oT inn^ T, inor. it i • u o stafs On 100-10 AU scalcs, using hnes of the CS molecule 

e.g., Blake et al. 1987, 1996; Turner 1989; Helmich & van , , , » , . , , ■ , . 

and continuum data. Adopting the resulting density pronles 



Dishoeck 1997; Gibb et al. 2000). This molecular complexity 
is thought to be associated with a specific stage in the evolu- 
tion of a massive young star — the so-called 'hot core' phase — 
in which selected molecules have just evaporated off the grains 



as input. Doty et al. (2002) calculated detailed gas temper- 
ature and chemical models to test these physical structures 
and, in particular, to study the chemical evolution of the high- 
mass star forming region AFGL 2591. Self-consistent models 
Send offprint requests to: P Stauber, e-mail: of the gas and dust thermal balance, chemistry, and radiative 
pascalst@astro.phys.ethz.ch transfer within YSO envelopes have also been performed by 
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Doty & Neufeld (1997) for high-mass objects and CeccarelH 
et al. (1996) for low-mass sources. Most recently, Rodgers & 
Charnley (2003) have modeled the chemistry of gravitationally 
collapsing low-mass YSO envelopes considering spatial den- 
sity variations and gradual heating of the gas as well as tempo- 
ral chemical evolution. 

All of these models do not consider a central source of 
UV radiation, and include (at most) an external radiation field 
comparable to the average interstellar radiation field. The only 
exception is the work by Ceccarelli et al. (1991), who mod- 
eled the chemistry around a low-mass YSO with central UV 
radiation and calculated radial abundances between w 1000- 
10^ AU. Their assumed density profile was flat rather than a 
power-law, whereby UV photons can penetrate further. The re- 
sults from such models resemble those from standard models of 
photon-dominated regions (PDRs), in which the UV radiation 
from a hot star impacts a nearby molecular cloud. In these re- 
gions the physical and chemical structure is determined by the 
incident FUV field (6 eV < /zv < 13.6 eV; see e.g., Sternberg & 
Dalgarno 1995, Hollenbach & Tielens 1999). PDRs are char- 
acterized by bright atomic fine-structure lines, strong high-/ 
CO emission, PAHs and specific molecules like CN, C2H or 
CO+ (e.g., Jansen et al. 1995, Storzer et al. 1996). Thus, while 
there is no doubt that UV photons affect their surroundings in 
the later stages of massive star evolution when the envelope 
is being dispersed, their importance in the early stages is still 
unclear. 

The question whether or not UV radiation plays a role in 
the chemistry of YSO envelopes has been a hotly debated sub- 
ject for some time. Because of the limited penetration depth of 
UV photons (typically a few magnitudes of visual extinction 
Ay), file bulk of the envelope with Ay > 100 mag must be un- 
affected. However, the presence of certain solid-state features, 
in particular the XCN or OCN" band at 4.62 yum, is often cited 
as an indicator of 'energetic processing' (e.g., Pendleton et al. 
1999). Whether this indeed implies the presence of UV radi- 
ation or simply requires higher temperatures is still uncertain. 
Indeed, recent laboratory experiments by van Broekhuizen et 
al. (2004) indicate that similar production efficiencies of OCN" 
can be obtained with thermal heating of HNCO-containing 
mixtures, reducing its diagnostic value. PAH emission is not 
detected in the earliest embedded stages, but this could be due 
either to freezing out of the PAHs onto grains or to the high ex- 
tinction at infrared wavelengths. The aim of this paper is there- 
fore to find other diagnostics of the inner UV field, in partic- 
ular from molecular tracers whose (sub-)millimeter emission 
can penetrate the entire envelope. 

We study here the chemistry in YSO envelopes under the 
influence of UV radiation from the central source as well as 
from an outside interstellar FUV radiation field. For this pur- 
pose, we have extended the time- and position-dependent de- 
tailed chemical model of Doty et al. (2002) to allow the impact 
of a central UV field on the inner envelope. In addition, we 
have used the resulting abundance profiles to compute emis- 
sion lines of selected species using the Monte Carlo code of 
Hogerheijde & van der Tak (2000), for direct comparison with 
observations. Although our applications are focused on high- 
mass YSOs, the results should be equally applicable to the en- 
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Fig. 1. Gas temperature distribution in the inner 200-250 AU for dif- 
ferent UV fluxes Go,in- The dust temperature r^,,,, equals the gas tem- 
perature for low UV fields (Go < 100). 



velopes around low-mass sources, albeit that the origin of the 
UV radiation may be slightly different. 

This paper is organized as follows. In Sect. 2, the physical 
and chemical model is described. The model is then applied 
to the envelope around AFGL 2591 and the results are dis- 
cussed. We examine the highly abundant molecules CO, H2O 
and CO2 in Sect. 3.1. In Sect. 3.2 we discuss the 'PDR-related' 
elements and ions C, O and C^. The reactive ions CO^, re- 
cently observed toward the high-mass star forming region W3 
IRS5 by Stauber et al. (2004), SO+, HCO+ and HOC+, which 
were also studied by Fuente et al. (2003) in PDRs, are exam- 
ined in Sect. 3.3. Hydrides and ions like OH, OH^, CH and 
CH^ are discussed in Sect. 3.4. These are important species for 
future space-borne observatories like Herschel where they will 
be observable with the high resolution spectrometer HIFI (de 
Graauw & Helmich 2001). In Sect. 3.5, CN -bearing molecules 
are studied since they are thought to be sensitive to UV radia- 
tion and readily observable with existing instrumentation. Line 
emission of CN and HC'^N is calculated and the results are 
presented and discussed in Sect. 4. We conclude by summariz- 
ing our results in Sect. 5. 



2. Model 

In order to calculate the abundances of each species at a certain 
time and distance from the central object, a physical and ther- 
mal model of the envelope is needed and initial conditions, i.e., 
initial abundances, have to be assumed. We adopt the values 
provided by Doty et al. (2002) and refer to that paper and ref- 
erences therein for more detailed information. However, a brief 
outline is given in this section. 
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' ' Following the arguments summarized in van der Tak et al. 

(1999), we assume a distance to AFGL 2591 of 1 kpc. 

2.2. Chemical model 

The chemical model is based upon the UMIST gas-phase 
chemical reaction network (Millar et al. 1997) and basically 
solves the equations for molecular evolution given by 
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Fig. 2. Column densities of various species as functions of Go.in! nor- 
malized to the value for Go.in = 0. 

2.1. Physical and thermal model 

Van der Tak et al. (1999, 2000) fitted their observations of 
AFGL 2591 with a spherical power-law density distribution 
n = miralrY («o = 5.8x 10'*cm-^ r,, = 2.7x 10"* AU, a = 1.0) 
by modeling CS emission lines over a large range of critical 
densities and fitting the continuum emission. Our model covers 
a region that extends from (rin, Wm) ~ (200 AU, 5.9 x 10*" cm"^) 
to (rout, "out) ~ (29 000 AU, 4.2 x lO^^cm^^). The dust tem- 
perature profile was determined from a self-consistent solu- 
tion of the continuum radiative transfer problem. As shown 
by Doty & Neufeld (1997) and Doty et al. (2002), the as- 
sumption that the gas and dust temperatures are equal, Tgas ~ 
Tdast, holds throughout most of the envelope because of the 
high density. However, at the inner edge, photoelectric heat- 
ing and other heating processes proportional to the radiation 
field (e.g., collisional de-excitation of vibrationally excited H2 
pumped by FUV radiation) can become important (Tielens & 
Hollenbach 1985, Sternberg & Dalgarno 1989). Comparison of 
these heating rates with the cooling rate due to gas-dust colli- 
sions (Hollenbach & McKee 1989) shows that the difiference in 
gas and dust temperatures is only a few K for low FUV fields 
(Go,in ^ 100; see Sect. 2.2 and 2.3 for definition). For higher 
FUV fields, the gas temperature becomes significantly higher, 
with values up to a few thousand K at the edge. 

Figure 1 illustrates the gas temperature in the inner re- 
gion computed taking the above-mentioned processes into ac- 
count for values of Go,in up to 10^. The gas temperature drops 
quickly within Ay < I and equals the dust temperature further 
in the cloud. The dust temperature — in particular the freeze- 
out radius (Fdust ~ 100 K) — is not affected beyond Ay x 1, 
since the FUV photons are quickly re-distributed in frequency 
and the total luminosity of the system is unchanged (Doty & 
Leung 1994). The physical structure is kept constant with time. 



dn(i) v-i v-i 



(1) 



j,k 



where n{i) is the number density of species /; kij and lijk 
are the rate coefficients for production and destruction of the 
given species. The main inputs of the model are the temper- 
ature and density distribution, the cosmic ray ionization rate 
f » 5.6 X 10 '^ s' for AFGL 2591, and the strength of the 
FUV field Go with respect to the average interstellar flux at 
6 < hv < 13.6 eV of 1.6 X lO^^ergcm^^s '. The model then 
calculates the time- and depth-dependent chemistry for 395 
species and 3865 reactions. The chemical abundances are cal- 
culated over a range of 34 radial grid points, providing a time- 
and space-dependent chemical evolution. The approximate ef- 
fects of freeze-out onto dust grains are included by initially 
depleting certain species below lOOK. The effects of changing 
the cosmic -ray ionization rate have been investigated by Doty 
et al. (2002) and are not further discussed in this paper. 

2.3. UV radiation field 

In our spherical model, the envelope can be irradiated by UV 
photons from the inside by a central source and from the out- 
side by the interstellar radiation field (ISRF). The UV flux is 
characterized in units of Go (see Sect. 2.2). We therefore denote 
the enhancement factor due to the inner source Go.in and that 
connected to the ISRF Go,out, assuming the same spectral en- 
ergy distribution. The rates for the photodissociation and pho- 
toionization reactions can then be approximated by (see also 
van Dishoeck 1988) 



k = Go.inCe '""'" + Go nutCe" 



yT„, 



l»-'l 



(2) 



where C and y are fit parameters and Tin and Tout are the at- 
tenuation factors, approximated by t « Ay, the total visual 
extinction Ay. The values of A^ at each depth are calculated 
from the total hydrogen column using the conversion factor 
A^(Htot) = 2A^(H2) + N(Ui) ~ 1.86 X lO^'Aycm-^, con- 
strained from observations of translucent clouds by Rachford 
et al. (2002). For simplicity we start our model with an atten- 
uation of Tin = at 200 AU and vary the inner UV flux from 
Go,in = to Go,in = 10^. We have also carried out different 
models by varying Go,out from up to 100 but the changes of 
the total radial column densities due to different ISRF fluxes are 
very small. The standard value for the ISRF flux in our models 
is therefore Go,out - 1 ■ 

To obtain a rough indication of the inner UV field, we as- 
sume a luminosity of 2 x 10^ L© and an effective temperature 
of 3 X 10"^ K for the central source as suggested by van der Tak 
et al. (1999). A large fraction of this luminosity may be due to 
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accretion (Osorio et al. 1999). This luminosity and stellar tem- 
perature would yield Go,in ~ 2x 10*^ at 200 AU in the absence of 
any absorption. Hence an estimated UV field strength of Go ~ 
10 corresponds to a dust optical depth t ~ 17 at this distance. A 
first question is whether such a relatively low value for t corre- 
sponds to a mass that is high enough to provide a typical mass 
accretion rate for high mass YSOs of ^ lO^'^-lO^^Moyr"' 
(e.g., Behrend & Maeder 2001, McKee & Tan 2003). Using 
the relation between the total visual extinction and the hydro- 
gen column density given earlier in this Sect., Ay ~ 17 mag 
corresponds to a column density of A^(H2) ~ 2 x 10^^ cm"^ 
or a density of n a; 10^ cm"^ at 200 AU. Assuming spherical 
free fall inflow and a central mass of M, = IOM0, this den- 
sity would be consistent with an accretion rate of (e.g., Shu 
et al. 1987) M,ec = 47rpP'\2GM,y'^ ^ IQ-^MQyv-K It is 
expected that these high accretion rates 'choke off' an incipi- 
ent H II region (Walmsley 1995) and that the hydrogen-ionizing 
photons (E > 13.6 eV) cannot travel far from the protostar be- 
fore being absorbed by infalling matter. The result would be 
a small Stromgren sphere around the central protostar, with 
the non-ionizing photons allowed to escape to larger distances. 
Both the luminosity and effective temperature of the radiation 
are actually expected to evolve during the protostellar stage 
(e.g., McKee & Tan 2003), but such detailed modeling is be- 
yond the scope of this study. Note also that the geometry of this 
dust and inflowing material is not known and is therefore not 
included in our envelope model. If our model were extended in- 
ward to account for this extra material, only the first few radial 
positions in our grid up to Ay ~ I would change. 

Our models are spherically symmetric and do not include 
geometrical effects like outflow cones or clumps in the central 
regions that could allow UV photons to travel further into the 
envelope. The effects of UV radiation escaping through outflow 
cones and scattering back on the surrounding envelope have 
been studied by Spaans et al. (1995) for the case of low-mass 
YSOs. Finally, it should be mentioned that the photoreaction 
rates (2) are fitted for plane parallel models but are applied to 
spherical models. In spherical symmetry, the FUV photons can 
penetrate further on average, but the gas temperature is some- 
what lower assuming the same density structure (Storzer et al. 
1996). The authors conclude that spherical geometry only be- 
comes important for clouds with n ^ 10^ cm ^^ exposed to 
FUV fields with intensities Go ~ 10^ and a H2 column density 
of less than » 10^^ cm"^. We therefore neglect these effects. 

3. Chemistry results 

We have modeled the envelope around AFGL 2591 with dif- 
ferent UV fluxes Go,in at the inner boundary and compared our 
results to those of Doty et al. (2002), representing the model 
Go,in = in this paper Photodissociation and photoionization 
processes are generally fast for the main species in the regions 
where they are significant (a; 10"*^ s"') so that the effects are 
not significantly different for 3x10-' and 3x10^ years. The re- 
sults are therefore presented for f = 3 x 10"* years, the chemical 
age of AFGL 2591 as proposed by Doty et al. (2002). Total col- 
umn densities for a selection of species are predicted in Table 
1 . The selection criteria are: (i) the species have a column den- 



sity larger than A^(X) = 10'" cm"^; (ii) the relative difference 
of the total column densities with and without an inner UV 
source is more than x 30%; and (iii) the species are observable 
with present or future instrumentation. In addition we show 
the abundances of the model predictions for basic molecules 
(other than H2) such as CO, H2O, CO2 and N2 at the bottom 
of Table 1 . In general, the species whose abundances are par- 
ticularly enhanced are the same as those found in traditional 
PDR models, i.e., radicals and ions, and their chemically re- 
lated species. The reason that most ratios given in Table 1 are 
> 1 is that small differences in very abundant species like CO or 
H2O can make large differences in much less abundant species. 

In Fig. 2, the column densities for various molecules as 
functions of Go,in are presented, normalized to the values for 
Go,in = 0. Figs. 3 to 6 show the depth-dependent fractional 
abundances for selected species in the inner 2500 AU. The ra- 
dial scale is taken to be somewhat smaller in these plots than the 
full extent of the envelope (29 000 AU, log(r) = 17.6) in order 
to emphasize the region which is influenced most by the inner 
UV field (see Doty et al. 2002 for variations in the outermost 
regions). Even for the highest values of Go,in, the effects of the 
UV radiation do not penetrate beyond the freeze-out radius at 
100K(log(r) = 16.2). 

Several molecules shown in Figs. 3 to 6 have local minima 
or maxima at certain positions (for example H2O at log(r) ss 
15.5 for Go.in = 10 or HCO+ at log(r) a 15.7 for Go.in = lO^). 
These peaks are not caused by numerical artifacts but are due 
to different chemical affects, some of which are explained in 
the following sections. 

Table 1 lists the radial column densities A^radiai = / "(O dr. 
For comparison of observed emission lines to model-predicted 
column densities, beam-averaged column densities defined by 
-^beam = ff n(z, p) G(p) IjT p dp dz / J G(p) IjT p dp havc to 
be calculated, where p is the impact parameter and G(p) is 
the beam response function. Beam averages are centered on 
the source whereas the calculated radial column densities are 
through the cloud center. A^beam, generally has slightly lower 
values than A^radiai for typical observing beams. Rather than 
calculating A^beam, however, we have computed the actual line 
fluxes by solving the equations of molecular excitation and 
radiative transfer throughout the envelope and convolving the 
emerging emission with the telescope beam (see Sect. 4). 

3.1. Basic molecules: CO, H2O, and CO2 

As seen in Table 1 and Fig. 3, the total radial column den- 
sities of the high abundance molecules CO, H2O and CO2 
are not affected significantly by an inner UV field. The major 
destruction mechanism for these molecules is photodissocia- 
tion. Measurable changes can only be seen for relatively high 
U V fields in the dense region within » 500 AU from the cen- 
tral source. For H2O another efficient destructive reaction is 
the production of HCO^ and HOC^ through reactions with 
C^. The slight increase in the relative abundance of H2O at 
r a 4x 10'^ cm for Go,in - 10 and Go,in = 100 is caused by the 
protonation of HCN through H3O+ which leads to HCNH+ 
and H2O, and by the high-temperature reaction of molecu- 
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Table 1. Predicted total radial column densities A'radiai for different UV fluxes Go.in- 



Species 


N [cm-- 


] Ratio 


N [cm--^] 


Ratio 


N [cm--^] 


Ratio 


N [cm-^] 


Ratio 




Go.in = 1 





Go,in = 100 




Go.i„ = 10' 




Go,n = 10^ 




H 


4.5(19) 


5.9 


1.3(20) 


17.6 


2.3(20) 


30.1 


1.4(21) 


187.1 


C 


1.0(17) 


95.2 


1.9(17) 


172.7 


4.4(17) 


400.8 


8.0(17) 


727.2 


N 


7.4(16) 


2.1 


1.4(17) 


3.9 


2.8(17) 


8.0 


4.9(17) 


14.0 





5.7(18) 


1.0 


6.1(18) 


1.1 


7.0(18) 


1.2 


8.8(18) 


1.6 


CH 


1.6(13) 


1.7 


3.7(13) 


3.9 


7.7(13) 


7.8 


1.7(14) 


18.2 


CH2 


1.2(14) 


59.3 


2.0(13) 


9.6 


1.7(13) 


8.3 


6.6(13) 


32.1 


CH3 


1.6(14) 


1.1 


1.6(14) 


1.1 


1.4(14) 


1.0 


1.4(14) 


1.0 


OH 


1.2(16) 


1.3 


1.6(16) 


1.7 


2.0(16) 


2.1 


1.2(16) 


1.3 


C2 


1.2(14) 


4.8 


1.0(14) 


3.9 


1.5(14) 


5.7 


1.7(14) 


6.8 


HCN 


1.4(17) 


238.6 


6.3(16) 


105.5 


8.2(15) 


13.7 


4.7(14) 


0.8 


HCO 


1.8(12) 


1.0 


1.7(12) 


0.9 


1.5(12) 


0.9 


1.3(12) 


0.7 


OCS 


4.0(14) 


1.5 


1.1(14) 


0.4 


7.1(13) 


0.3 


4.9(13) 


0.2 


CN 


3.1(14) 


5.6 


4.2(14) 


7.5 


1.8(14) 


3.3 


9.2(13) 


1.7 


cs 


1.0(15) 


2.4 


7.9(14) 


1.9 


4.2(14) 


1.0 


3.7(14) 


0.9 


HCS 


3.3(14) 


34.1 


1.1(13) 


1.1 


1.4(11) 


l(-2) 


1.5(11) 


2(-2) 


C2H 


1.5(15) 


221.4 


1.4(14) 


20.5 


8.5(13) 


12.7 


7.3(13) 


10.8 


C3 


2.9(13) 


1.6 


3.0(13) 


1.7 


2.9(13) 


1.7 


2.2(13) 


1.2 


C3N 


2.3(13) 


1.4 


1.8(13) 


1.1 


1.7(13) 


1.0 


1.6(13) 


1.0 


C5 


7.0(12) 


1.3 


7.3(12) 


1.4 


6.8(12) 


1.3 


5.9(12) 


1.1 


C5N 


6.3(11) 


2.0 


3.2(11) 


1.1 


2.2(11) 


0.7 


1.9(11) 


0.6 


HNC 


2.0(15) 


2.7 


1.2(15) 


1.6 


7.3(14) 


1.0 


6.5(14) 


0.9 


HC3N 


4.6(14) 


1.6 


3.0(14) 


1.0 


2.9(14) 


1.0 


2.9(14) 


1.0 


C2S 


6.6(14) 


14.0 


1.5(14) 


3.2 


9.6(13) 


2.1 


1.2(14) 


2.6 


C3S 


2.6(13) 


4.2 


1.5(13) 


2.4 


1.1(13) 


1.8 


6.6(12) 


1.1 


CH2CN 


1.3(16) 


1074.3 


6.8(15) 


566.0 


6.2(14) 


51.9 


7.4(12) 


0.6 


CH3CN 


8.9(13) 
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a(b) means a x 10*. The ratio is the enhancement relative to the model without an inner UV source, i.e. A'(Go,in = 10^)/A'(Go,in = 0) 



lar hydrogen with OH, also producing H2O. All reactants are 
slightly increased by a moderate UV field, enough to speed 
up these reactions. At x; 5 x 10'^ cm the relative abundance 
of CO is the same for all UV fields. This is slightly closer to 
the star than the enhancement of OH for high UV fields (see 
also Fig. 5) and since CO + OH — > CO2 + H is the dominant 
reaction to produce gas-phase CO2, the enhancement of OH 
leads to the small increase of CO2 between 5-6 x 10'^ cm for 



Go.in ~ 10^. The sudden decrease of the H2O and CO2 abun- 
dances at w 2 X let'em (r w 100 K) is due to the effects of 
their freeze-out onto grains (see Doty et al. 2002 for more in- 
formation). The abundance minima of H2O, CO and CO2 for 
Go,in > 10^ are due to the higher gas temperature for these high 
FUV fields: the abundances are rising again toward the proto- 
star. 
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Fig. 3. Depth-dependent fractional abundances of the basic molecules CO, H2O and CO2 and the 'PDR-related' species O, C and C^ for 
different UV fluxes. 



CO2 is a highly abundant molecule in our model with a col- 
umn density of the order of lO'^ cm"^. Boonman et al. (2003) 
derive a total column density of 2.6 x 10"" cm"^ for AFGL 259 1 
from observations with the Infrared Space Observatory. Even 
high UV fields do not reduce much CO2 . This overprediction 
of CO2 has also been discussed by Doty et al. (2002), and our 
new models add little to this discussion: enhanced UV pho- 
todissociation does not solve the problem. 



3.2. PDR-related species: O, C, and C+ 

Fig. 3 also shows the dominant line emitters, and thus coolants, 
in a PDR: O, C and C+. Photodissociation of CO is the 
main source for both carbon and oxygen. The major destruc- 



tion mechanism for carbon in a low (Go,in = 10) UV field is 
the reaction with OH leading to CO and atomic hydrogen. 
Photoionization of carbon becomes important for higher UV 
fields and is the main production channel of C^. The some- 
what unexpected jump of C^ at * 2 x 10'* cm can be explained 
by the freeze-out of H2O: the main destroyer of C^ at this dis- 
tance is depleted onto grains below lOOK. The relatively high 
atomic oxygen abundance for T < 100 K is due to our initial 
elemental abundance, consistent with Meyer et al. (1998) for 
diffuse clouds. 

The main production channel for O^ is not the photoioniza- 
tion of atomic oxygen (which cannot occur for photon energies 
less than 13.6 eV) but the reactions H+ -H O ^ 0+ -H H and C+ 

-H O2 ^ 0+ H- CO. 
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The observed column density of atomic carbon toward 
AFGL 2591 (A^ < 6.8 x lO'^ cm"^ see Doty et al. 2002) can be 
reproduced with a UV field of Go,in > 10. The high C^ abun- 
dances (A^ < 6.8 X lO'^ cm"^) are consistent with a UV field of 
Go,in > 10^. However, as discussed in Sect. 4, the observed an- 
tenna temperatures require high abundances of C and C^ to be 
extended over a larger region than just the inner few hundred 
AU. Thus, these species likely trace the outer radiation field, 
Go.out ~ 10 - 100, rather than an inner UV field. 



3.3. Reactive ions 

From Table 1, it is seen that HOC^ is much more aff'ected by 
the strength of the UV field than HCO^ . The reason for this 
is that in our models HOC^ is basically formed by the two re- 
actions H+ -H CO ^ HOC+ + H2 and C+ + H2O -^ HOC+ 
+ H, which are both highly dependent on the ionization frac- 
tion and thus on the incoming UV field. In contrast more pro- 
duction channels are possible for HCO^, e.g., HOC^ + H2 
-^ HCO+ + H2, thus more HCO+ than HOC+ is produced. 
Another reason for the much higher abundances of HCO^ is 
the jump at r = 100 K where H2O freezes out: since HCO^ 
+ H2O — > H30^ + CO is the main destruction channel, much 
more HCO^ remains in the region below lOOK. Water is im- 
portant, however, to produce HOC^, hence the HOC^ abun- 
dance is lower once water has depleted onto grains. 

The HCO^/HOC^ abundance ratios predicted by our 
chemical models are in the range of a; 10^-10^ which is at 



the upper limit of the ratios observed toward PDRs. Fuente et 
al. (2003) reported ratios of 50-120 toward the reflection neb- 
ula NGC 7023 — the lowest ratio measured so far — up to 
450 toward the planetary nebula NGC 7027. Our models only 
reproduce such low values for Go,in ^ 10^. However, the ob- 
served peak values of HOC^ toward several PDRs are in the 
range of the predicted column densities for Go,in - 10-100. 

Van der Tak et al. (1999) find HCO+ abundances of the 
order of w 10"^, which are well fitted by our models with 
Go.in = 10. Unfortunately, there are no observations of HOC^ 
toward AFGL 2591 to date. The question therefore is whether 
these low ratios also hold in massive star-forming regions. 

For low UV fields the reactive ion CO^ is produced mainly 
by He+ + CO2 ^ CO+ H- O -H He and C+ + CO2 ^ CO+ 
+ CO. For higher UV fields the reaction C+ H- OH -» CO+ + 
H is the most efficient to build up CO^ . Its destruction for all 
UV fields is mainly by reactions with H2 . The most efficient 
pathway to SO+ is S+ + OH -^ SO++ H. Its destruction is 
caused mainly by dissociative recombinations with electrons. 



3.4. Hydrides 

As in traditional PDR models, the abundances of unsaturated 
hydrides are strongly increased by UV radiation. OH and CH3 
are the species that are the least affected by the inner UV 
field. CH^ on the other hand is enhanced by more than a fac- 
tor of 10^ even for a moderate UV field and its column den- 
sity grows steadily with increasing UV field strength. To date. 
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there are observations only of OH for AFGL 2591 reporting 
A^ > 4.7 X 10"* (see Doty et al. 2002). This column density is 
well fitted by our models for the region r < 100 K. Our models, 
however, produce much more OH in the warmer, inner region 
(see Fig. 5). Without UV radiation, OH is mainly produced 
through collisions of atomic oxygen with molecular hydrogen 
in the inner warm part. In the presence of a UV field, photodis- 
sociation of H2O becomes the dominant production route. The 
OH destruction is primarily through C^ , and for very high UV 
fields (Go,in - 10^) the photodissociation of OH becomes im- 
portant too. For UV fields with Go,in > 10^, CH^ is efficiently 
produced in the warm inner regions by reactions of C^ with 
molecular hydrogen. For smaller UV fields, CH^ is produced 
mainly by the reaction HCO^ + C ^ CH+ + CO. For higher 
UV fields, photoionization of CH becomes the main mecha- 
nism for CH^ production. In all our models, CH^ is mainly 
destroyed in collisions with H2. 

For smaller UV fields, CH is mainly destroyed through CH 
H- H2 — > CH2 + H and built up by the reverse reaction. OH^ 
is primarily produced by collisions of singly ionized oxygen 
with molecular hydrogen. The jump at » 2 x lO'*" cm is due to 
the increase of atomic oxygen. Charge exchange reactions then 
enlarge the OH^ abundances, in particular the reaction H3 + 
O — > OH^ -I- H2 . OH^ is most efficiently reduced in collisions 
with H2. OH and OH^ play a crucial role in the gas-phase 
chemistry of water. 

Several species such as C^ and CH^ are not only excellent 
tracers of enhanced UV fields but also of the strength Go since 
their abundances rise steadily with increasing field strengths. 
The fact that C^ can also be enhanced by an outer UV field re- 
duces its value as a diagnostic, however. An alternative tracer 
is the CH^/CH ratio, which is 10^-10^ times higher if the 
inner UV is included. The column density ratio is shown in 
Fig. 7 as a function of Go,in- Similar arguments hold for other 



hydrides such as CH2 and CHj which are less easily observ- 
able. Hydrides like CH and CH^ are also important for building 
complex hydrocarbons (see e.g., van Dishoeck & Hogerheijde 
1999); indeed, a rich carbon chemistry has been observed in 
PDR layers (e.g., Hogerheijde et al. 1995, Teyssier et al. 2004). 

It is somewhat surprising that CH3 is much less affected 
by the inner UV field than, e.g., CH2. This is only partly true 
since the production and destruction mechanisms of CH3 do 
change by the impact of the inner UV field. Without the inner 
UV source CH3 is efficiently built up by the reaction He^ -H 
CH3OH -^ OH+ + CH3 + He. With the onset of an inner UV 
field, CH3 is mainly the photodissociation product of species 
like CH3CN, CH4 and CH3OH. CH3 forms H2CO when it re- 
acts with O, which is the main destruction channel. Since an 
inner UV field clearly enhances atomic oxygen (see Fig. 3), 
this reaction is very efficient deeper in the cloud where CH3 is 
destroyed fast. Thus, although CH3 is first enhanced in the very 
inner region, it becomes less abundant deeper in the cloud. This 
leads to an average radial column density more or less equal to 
that where no inner UV radiation is considered. 



3.5. CN -bearing species 

The chemistry of CN and HCN is thought to be strongly af- 
fected by UV radiation and the CN/HCN ratio has been pro- 
posed to be a good tracer of PDR chemistry (Fuente et al. 
1995, Jansen et al. 1995). Since our envelopes are more ex- 
tended than normal PDRs, we predict lower CN/HCN ratios 
than those observed in PDRs where the ratio is typically ss 1 
in the shielded regions and > 1 in regions exposed to UV 
radiation. Observations of CN and HCN toward AFGL 2591 
confirm these low ratios. The observed abundance of CN (see 
Doty et al. 2002) is well fitted with Go.in = 10-100. The HCN 
chemistry is discussed in detail by Doty et al. (2002), espe- 
cially the situation in the warm gas-phase (T > 230 K) region 
where a jump has been suggested from x(HCN) - 10"^ to 
x(HCN) = 10^^ by Boonman et al. (2001). The observed HCN 
column density is reproduced in our models for Go,in = 10 
within a factor of 2, whereas it is underproduced in models 
without an inner UV field. Alternatively, the model could be ex- 
tended inward to the radius where high temperature chemistry 
may enhance the HCN abundances. However, a closer position 
to the star also means an increased UV field, and HCN is more 
vulnerable to photodissociation. UV radiation may therefore 
be an explanation for the relatively high observed HCN abun- 
dances. The CN/HCN ratio is only enhanced for Go,in ~ 10^. 
For lower UV fields, the ratio is even smaller than that without 
UV irradiation (see also Fig. 7). Hence, the CN/HCN ratio is 
not a good UV tracer for this high mass star-forming region, in 
contrast with the situation for PDRs. 

Fig. 6 shows the radial dependences of the abundances of 
CN, HCN, HNC and the reactive molecular ion HCNH+. In 
the warm high density regime, the production of HCN through 
collisions of CN with H2 is slightly more efficient than pho- 
todissociation of HCN, hence more HCN is produced than CN. 
Collisions of HCN with H30^, however, produce water and 
HCNH^ , which can recombine to form HNC or CN . Since this 
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is the main path to HNC in our models, more HCN is produced 
compared to HNC. CH3CN is generated through H4C2N+, 
whereas the association of CN with CH3 plays a minor role. 
CH3CN and HNC are also photodissociated to CN. No grain 
surface production and subsequent evaporation of CH3CN are 
assumed. The following network is valid for UV-affected CN 
chemistry in our models: 



HNC/HCN= 

field. 



0.01, which is in good agreement with a low UV 
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The HNC/HCN ratio is * 0.01-0.03 for low UV fields. The 
ratio is approximately unity for Go,in = and slightly higher 
(1.4) for high UV fields (Go.in = 10^). The variation of the 
ratio is mainly due to the changing HCN abundances. The ob- 
served HNC/HCN ratio varies strongly from source to source 
but is generally < 1. For AFGL 2591, the observed ratio is 



4. Calculated line emission 

The temperature, density and abundance profiles presented in 
Sect. 3 have been used as input to the excitation and Monte 
Carlo radiative transfer code of Hogerheijde & van der Tak 
(2000) to compute the emerging line intensities. The intrinsic 
(turbulent) line profile is taken to be a Gaussian with a Doppler 
parameter (1 /e width) of 1 .6 km s"' , independent of radius. The 
line profiles are convolved to an appropriate telescope beam 
size. The low-lying transitions trace primarily the low density 
outer YSO envelope whereas the higher-lying lines are only ex- 
cited in the high density inner envelope. Thus, the intensities of 
the higher-lying transitions of species like CN or HCN, whose 
abundances are affected by UV irradiation, are expected to be 
enhanced and the effect on the lines ratios should be observ- 
able. 
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4.1. CN and HCN emission lines 

Fig. 8 shows emission lines of CN for different inner UV fields. 
While the lowest transition CN 1-0 line has the same intensity 
for all UV fields, the higher transitions show clear differences 
for different UV strengths. The CN 3-2 line intensity is in- 
creased by a factor of x; 2 for Go,in = 10 and Go,in = 100. The 
enhancement is even more noticeable in the CN 4-3 transition. 
The lines for Go,in = 10 and Go,in = 100, respectively, are a: 3 
times stronger than that for Go,in = 10^, which in turn is approx- 
imately twice as strong as that for Go,in = 10^ and Go,in = 0. 
Since these enhancements are clearly larger than the typical 
30% observational errors, the impact of central UV photons 
on the inner envelope is predicted to have measurable conse- 
quences for higher CN transitions. The line ratios for Go,in = 
are CN 1-0/3-2/4-3 ~ 100/11/2, whereas for Go.in = 10 they 
are approximately 100/18/12 and for Go,™ = 100 CN 1-0/3- 
2/4-3-100/17/11. 

In Fig. 9, it is seen that the higher-lying HC'^N lines are 
similarly enhanced for a moderate UV field with Go,in = 10. 
The optically thin HC'^N isotope is modeled, since the syn- 
thetic HCN lines are optically thick and therefore not ideal to 
show the effects of inner UV enhancement (see also the dis- 
cussion by van der Tak et al. 1999 on HCN modeling). The 
HCN/HC'^N ratio is taken to be 450 (Wilson & Rood 1994). 
The 7 = 3-2 line is enhanced by a factor of 7 and the 7 = 4-3 
line a factor of w 26 for moderate UV fields. For the low- 
est transition, all line intensities are the same, except that for 
Go.in = 10, which is approximately 9% higher and thus not dis- 
tinguishable from the other lines if calibration errors are taken 
into account. The intensity ratios for Go,in = are HC'^N 1- 
0/3-2/4-3 ~ 10/10/4, for Go.in = 10 HC^^N 1-0/3-2/4-3 ^ 
10/64/103 and for Go.in = 100 they are approximately 10/41/70. 

Line intensities have also been calculated for an outer UV 
field of Go,in = 10-100. It is found, however, that the CN and 
HCN lines are not very sensitive to the outer UV field due to 
the much lower densities in those regions, with enhancements 
in Une intensities of only 10-20%. Thus, the CN and HCN 
lines trace the central UV field rather than the outer interstellar 
radiation field. This is in contrast with species such as C and 
C^, whose transitions are readily excited at low densities in the 
outer envelope (see below). High spatial resolution interferom- 
etry observations of both low- and high-J CN and HCN lines 
should be able to directly constrain the extent over which the 
inner and outer UV field affect the chemistry. 

4.2. C and C^ emission lines 

In addition to the molecular lines mentioned above, the line 
profiles of the fine-structure transitions of C and C^ have 
been calculated. Although the abundances of both species are 
strongly enhanced by the inner UV field, the resulting antenna 
temperatures are much lower than observed by van der Tak et 
al. (1999). Thus, even though the column densities including 
an inner UV field match those derived from observations, the 
beam dilution of the inner region is so large that the result- 
ing fluxes fall short by 1-2 orders of magnitude. Therefore, 
models have also been run in which the outer radiation field is 
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enhanced above the average ISRF. Such outer enhancements 
could be due to other nearby young stars. For example, for 
AFGL 2591 a nearby Hii region (6000 AU projected distance) 
has been found (see e.g., van der Tak et al. 1999, Trinidad et al. 
2003), whose exciting star could illuminate (part of) the outer 
AFGL 2591 envelope. An outer UV field with Go,out = 100 
and an enhanced gas temperature due to photoelectric heat- 
ing matches the [C I] 492 GHz observations best. Since C^ 
is mainly produced through photoionization of atomic carbon 
(ionization energy 1 1 .26 e V), the same issue may be true for 
the [C II] fine. 

5. Summary and conclusions 

We have extended the detailed chemical models of Doty et al. 
(2002) to study the influence of UV irradiation from AFGL 
2591 on the chemistry in the surrounding envelope on 200- 
29 000 AU scales. From Figs. 3 to 6, it is seen that an inner UV 
flux affects only the region within 500-600 AU from the star. 
The species whose abundances are most enhanced are radicals 
and ions, similar to the case of normal PDRs. However, the 
chemistry differs from that in normal PDRs in various details 
owing to the higher temperature and higher H2O abundance in 
the inner YSO envelope. 

By comparing our results to those of Doty et al. (2002), 
we find that agreement with the observed column densities is 
improved with a moderate UV field of Go,in = 10-100. This 
is seen quantitatively in Fig. 10, which shows the absolute 
differences between the model results and the observations, 
A = F liZi I log(A^obs) - log(A^mod)l, summed over afl species. 
Species like CN and HCN are particularly improved with a 
modest UV field. However, the high temperature chemistry of 
HCN in the inner region still needs further investigation. The 
CN/HCN ratio is found not to be a good tracer for Go,m for 
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AFGL 2591. For low UV fields, the extra CN reacts with H2 
to HCN in the warm, dense region thus enhancing HCN more 
than CN, in contrast to typical PDRs. However, enhanced inten- 
sities of the 3-2 and 4-3 lines of CN and HC'^N are predicted 
compared with the lower 7=1-0 lines, which should be ob- 
servable both with single dishes and with submillimeter inter- 
ferometers. Another strong indicator for enhanced UV fields is 
the CH^/CH ratio with values between 10"^ and 10 '. Neither 
species can readily be observed with ground based observato- 
ries, but are excellent targets for Herschel-HIFI. 

Most other species have their maximum column density at 
lower UV fields (see Table 1). This is due to different preferred 
chemical networks for each strength of the UV field; in par- 
ticular, many species are photodissociated away for very high 
UV fields. Like in Doty et al. (2002), the CO2 column density 
is overpredicted by a factor of 40 in our model. In this sense, 
UV irradiation from the central star proves not to be an efficient 
destruction mechanism for CO2 . Impulsive heating events and 
X-ray chemistry may destroy carbon dioxide more efficiently. 

Our model results indicate strongly that UV radiation from 
the high mass YSO AFGL 2591 cannot be neglected in the sur- 
rounding envelope chemistry. The inferred value of Go,in = 10- 
100 for an inner UV field is reasonable for a young massive 
star like AFGL 2591 and the corresponding column density 
and dust opacity can be justified by the mass inflow rate even 
though we do not present a fully self-consistent model (see 
Sect. 2.3). Our models assume a spherical symmetry and ne- 
glect geometric effects. High spatial resolution data of species 
like CN should also be able to establish the importance of non- 
spherical effects, such as the escape of UV photons through the 
outflow cones. 

In recent years it has become clear that YSOs are 
also strong X-ray emitters. The observed X-ray luminosities 
range from approximately 10^*^'^ to lO-'-'ergs"' (Feigelson & 
Montmerle 1999, Preibisch 1998). X-rays can also penetrate 
deeper into clouds than UV photons due to smaller absorp- 
tion cross sections at high energies. Most recently Doty et al. 
(2004) suggested that X-rays may be responsible for a higher 
ionization rate in the low-mass hot core IRAS 16 293-2422. It 
is therefore important to also investigate the effects of X-rays 
on the chemistry in the envelopes of YSOs and such work is in 
progress (Stauber et al. 2004, in preparation). 
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